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We describe a new and consistent perturbation theory for solid-state 

quantum computation with many qubits. The errors in the implementation 
of simple quantum logic operations caused by non-resonant transitions are es- 
timated. We verify our perturbation approach using exact numerical solution 
for relatively small (L = 10) number of qubits. A preferred range of pa- 
rameters is found in which the errors in processing quantum information are 
small. Our results are needed for experimental testing of scalable solid-state 
quantum computers. 
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Several proposals for scalable solid-state quantum computers have been recently pub- 
lished [|ilHjT|]- For the most effective quantum information processing all of these proposals 
require operations with large number of qubits. In Ref. |E| quantum logic operations be- 
tween remote qubits were simulated in a chain of L = 1000 nuclear spins. Such a quantum 
problem can be solved only using approximate methods since the dimensionality of the 
Hilbert space increases as 2 L . The approximate procedure used in |12| was based on the 
selection of states generated as a result of resonant or near-resonant transitions, while the 
other transitions were neglected in a controlled way. In this Letter we develop a consistent 
approach to this problem based on perturbation theory. We use our procedure to analytically 
estimate the probability of generation of unwanted states caused by non-resonant transitions 
and verify that these agree with the exact numerical solution for the same problem with rel- 
atively small number of qubits (L = 10), for which the dimensionality of the Hilbert space 
is not very large (N = 2 10 = 1024). 
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Dynamics of a spin chain We consider a chain of identical nuclear spins placed in an 
external high magnetic field, B(x), which has uniform gradient in the direction of the chain, 
x. The nuclear magnetic resonance (NMR) frequency for the kth spin is u>k = 7-Bfc, where 
7 is the nuclear gyromagnetic ratio and B& is the ^-component of the magnetic field at the 
location of the kth spin. The gradient of the magnetic field provides a change in uj^ by the 



value 8uj between the neighboring spins. (For physical parameters see Ref. [12|.) 
The Hamiltonian of the spin chain in an external radio-frequency (rf) field is, 

H (n) = - 5>*jj - 2J J2 IZ Ai- 

k=0 k=0 

L-1 

e„(*)(n„/2) £[4~ex P (_^ nt + ^ + i+ e x P (i(u n t + <p n ))] = h + v {n \t), (i) 

k=0 

where J is the Ising interaction constant between neighboring spins and I| is the operator 
of the ^-component of spin 1/2; Q n , v n and tp n are the Rabi frequency, the frequency and 
the phase of the nth pulse, = 1% ± i^, and Q n (t) equals 1 only during the nth pulse. 

In this Letter we estimate the errors generated during the creation of the entangled 
state for remote qubits from the ground state of the chain by applying a single n/2 pulse 
and a sequence of it pulses in the system described by the Hamiltonian ([]]). First, we 
describe schematically the quantum protocol. The first n/2 pulse creates a superposition 
of two states with equal probabilities, |00 . . . 00) -> (l/ v / 2)(|00 . . . 00) + z|10 . . . 00)). The 
subsequent pulses, which we describe by the the unitary transformation, U, transform this 
state to an entangled state for remote qubits, 

*7-J=(|00...00) +z|10...00)) = -=(e < * 1 (|00...00) + e^ 2 |10 . . . 01)), (2) 
v 2 v 2 



where <j)\ and (j>2 are known phases (see Ref. |12|). The operator U realizes a particular case 



of the well-known CONTROL-NOT gate. It has the following properties 

-L|o 1 _ 1 o I _ 2 ...o 1 o ) = -L 
^IU- 1 o i _ 2 ...o 1 o„) = -L 



C/-^l°i-i i-! ■ ■ ■ °i°o> = 4?e'*' |0 L -,0i- 2 . . . 0,Oo>, (3) 



X 1 

t/-=|U_ 1 O i _ 2 ...OiOo> = - 7 =e<^-^)|l i _ 1 z .-2 • • -Oilo). (4) 



To accomplish the operation (|4]) in the system described by the Hamiltonian (|lj), we choose 
a sequence of 7r pulses with resonant frequencies which will be described elsewhere. If we 
apply the same protocol to the ground state (operation (|3|)), then with some probability 
the system will remain in this state because these pulses have the detunings from resonant 
transitions, A n 7^ 0. The near-resonant transitions have the probabilities |12| ], 

e n = (^ n /A n ) 2 sin 2 (A n r n /2), (5) 



where A n = y A n + Vt n is the effective field in frequency units and r n is the duration of the 
nth pulse. The values of the detunings are the same for all pulses, A n = A = 2J, except 
for the fourth pulse, where A4 = 4J. In our calculations we assumed Q n = Q for n 7^ 4 and 
f2 4 = 2Q, so that for all pulses the values of e n are the same, e n = e. 

We write the wave function, in the time-interval of the nth pulse, in the laboratory 
system of coordinates in the form, 



^(t) = exp 



L-l 



[u n t + <p n ) IZ k 

k=0 



% ot (t) = j2Mt)\p)eM-*x { ; ) t + & ) ), (6) 



where ^/ rot (t) is the wave function in a frame rotating with the frequency v n . xi = 
— (v n /2) Z)fc=o ^li ^ = J2k=o ^li a k = ~ 1 if t ne kth. spin of the state \p) is in the 
position 1 and a\ = 1 if the kth spin is in the position 0, \p) is the eigenfunction of the 
Hamiltonian H . The dynamics during the nth pulse with the frequency v n is described by 
the following Schrodinger equation for the coefficients A p (t), 

iA p {t) = (E p - Xp n) )A p (t) -§£ A p ,(t), (7) 

p' 

where the sum is taken over the states \p') connected by a single-spin transition with the state 
\p). Eq. © can be written in the form iA p (t) = H [ $A p ,{t), where H ( p ) = H$ - xJ, n) <W, 
where 5 PtP > is the Kronecker 5-symbol. 

Under the condition A <C 8uj the energy separation between the pth and the mth diagonal 

elements of the matrix H^l) connected by the resonant or near-resonant transition is — 

pp j 1 v 

£ffl\ = A^ m , where E^ 1 ' = E p — x p n ^- This is much less than the energy separation between 
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these diagonal elements and diagonal elements of the other states connected with the states 
\p) and \m) by the non-resonant transitions. In this case one can neglect the interaction of 
the pth state with all states except the state \m). In this approximation the Hamiltonian 
matrix TC p p, breaks up into N/2 approximately independent matrices 2x2, where N = 2 L . 

Errors in the creation of an entangled state for remote qubits. The matrix approach 
allows us to estimate errors in the logic gate @ caused by near-resonant and non-resonant 
transitions. Suppose that initially the eigenstate \p) of the Hamiltonian H is populated. 
We want to calculate the probability of non-transition to the state \r) with \E P — E r \ ~ Sou, 
where the states E p and E r are connected by a flip of k'th spin, whose NMR frequency 
differs by the value ~ Su from the frequency of the resonant transition. Since the matrix 
elements are small, we can write 

^) = < „, + £<_* ^W, ( 8 ) 

p Oq Op 

where prime in the sum means that the term with p = q is omitted; v)ff is the matrix 
element for transition between the states ^l"' and ipp^ n \ 

Because the matrix Ti^ is divided into 2 L_1 relatively independent 2x2 blocks, the 
wave function, ipp( n \ in Eq. (Q) is an eigenfunction of a single block 2x2 with all other 
elements being equal to zero. In the case for which the two states \p') and \p) are connected 
by a near-resonant transition, the function, ipp^ n \ in Eq. (§) has the form, ip®^ « [1 — 
(Q 2 /8A 2 )}\p) + (Q/2A)\p') . On the other hand, if the states \p) and \p') are connected by an 
exact resonant transition, we have = (l/y / 2)(|p) + \p')) and Vy = (l/\/2)(|p) — \p')). 
In both cases the probability of non- resonant transition from the state \q) to the state \p) 
connected with the state \q) by flip of the k'th spin (up to the value /i = (fl/u) 2 ) is 

where \k n — k'\ is the distance from the fc'th spin (whose NMR frequency, toy, is non-resonant) 
to the fc n th spin with resonant (or near- resonant) NMR frequency. 

The probability fiL-i (here the subscript of fi stands for the number of the resonant 
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spin) of generation of unwanted states by the first 7r/2 pulse in the result of non-resonant 
transitions is 



2 L-2 



1 



Ay E . 



(10) 



After the first ir/2 pulse, the probability of the correct procedure in implementation of the 
logic gate is Pi = 1 — fii-i- The probability of correct implementation of the operation @ 
by applying 2L — 2 pulses is 



1 



L-3 



P2L-2 = ^(1 - - - e)(l - 4/i L _ 2 - e)(l - Ho - e) J[ (1 - m - ef+ 

~ i=i 

-1 L-3 



(ii) 



8=0 



where the first term is provided by the operation (|3|), and the last term is due to 
(a) 
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FIG. 1. (a) The probability, T 3 , of generation of unwanted states in implementation of the logic 
gate (0) at e = 0. Filled circles connected by the solid line are the numerical results, dashed line is 
the analytic estimate; O = 2J/ v / (4k 2 - 1), k = 8. (b) P function of 0. Filled circles are the 
numerical results. The dashed line is the analytic estimate, 5uj = 100; J = 1, L = 10, (p n = for 
all n. 

Numerical results. In Fig. 1 (a) we compare the total probability, V = 1 — P2L-2, of 
generation of unwanted states with the result of exact numerical solution of Eq. (|7|) for 
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the case in which the probability of near-resonant transitions is negligibly small, i.e. when 
e = 0. From Fig. 1 (a) one can see that in this situation V decreases as 5ui increases. 

When e is large, e 3> Li, and A <C 5u, the probability of generation of unwanted states 
is mostly defined by the value of e and is almost independent of 8u. In Fig. 1 (b) we plot 
the probability V as a function of Q. The value of Su fixes the values of the minima, V m i n , 
in Fig. 1 (b): for larger 5u the minima in the plot in Fig. 1 (b) become deeper. The 
values of different minima in Fig. 1 (b) indicate the contribution of non-resonant processes 
to the probability V. Since the value of Q in Fig. 1 (b) does not change significantly, 
the contribution of non-resonant processes to the probability of errors is approximately the 
same for all Q and is equal to V m in- One can see that this contribution is negligibly small in 
comparison with the contribution of near-resonant processes (defined by e) for all Q, except 
for the small regions of Q, where e, defined by Eq. (|5|), is minimal. 
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FIG. 2. The total probability of unwanted states in implementation of the logic gate @ as 
a function of Su when the value of 8lo is comparable with the value of the detuning, A. (a) 
Q = 2J/V4/c 2 - 1, k=8 (e = 0), (b) Q, = 0.15 (e = 0.0039). Filled circles are the numerical results, 
solid line is the analytical estimate for V; J = 1, L = 10, ip n = for all n. 



When the perturbation parameters of the problem, e and /a, are small, using Eq. (|Tl|) one 
can easily estimate the probability of generation of unwanted states in the implementation 
of the logic gate (0) when the number of qubits in the spin chain is large (for example, for 



L ~ 1000). In this case the Eq. ( |TTD is important for estimation of errors, since the exact 
solution of the problem requires diagonalization of enormous matrices of size 2 L x 2 L . 

We should note that one more condition (except for e, \i <C 1) must be satisfied for 
Eq. (IT]) to be valid. The value of detuning, A, should be small in comparison with the 
difference between NMR frequencies of the spins, A <C 5uo. In Figs. 2 (a) and (b) we plot 
the probability, V, as a function of 8uj for small and large values of e. One can see that our 
results are valid only when A <C Su, in spite of the fact that the parameter A/5u does not 
appear explicitly in Eq. (|TT|) . From Fig. 2 (b) one can see that the probability of unwanted 
states, V for e ^> \i and for large 8uj becomes relatively independent of 8u. In this case the 
value of V is defined by e which, due to Eq. (||), does not depend on 8uj. 
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